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ABSTRACT 

Stability  in  the  sense  of  Lyapunov  is  investigated  for  systems 
of  ordinary  differential  equations.  Various  established  results 
are  discussed  and  a  new  set  of  sufficient  conditions  for  a  stability 
comparison  theorem  is  obtained.  The  treatment  is  principally 
concerned  with  linear  systems,  though  some  results  for  non-linear 
systems  are  considered  briefly. 
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I.   INTRODUCTION 

Stability  was  originally  a  physical  concept,  believed  first  to  have 
been  considered  by  Lagrange  in  connection  with  the  equilibrium  states 
of  conservative  systems.  Equilibrium  states  are  those  in  which  the 
velocities  and  accelerations  are  zero,  describinn  a  state  of  rest. 
A  system  in  equilibrium  may  react  in  one  of  two  manners  when  influ- 
enced externally;  it  may  exhibit  motions  which  remain  in  a  small 
reqion  about  the  equilibrium  position  (which  does  not  oreclude  the 
oossibility  that  the  system  return  to  its  state  of  equilibrium), 
or  it  may  exhibit  motions  which  cause  the  system  to  increasinqly 
deviate  from  its  equilibrium  state,  never  returninq  to  it.  In  the 
former  case,  the  equilibrium  state  is  said  to  be  stable,  while  the 
latter  case  describes  an  unstable  state  of  equilibrium. 

The  above  remarks  are,  of  course,  not  rigorous  and  are  intended 
solely  as  an  intuitive  guide  to  the  meaning  of  stability  of  equili- 
brium states  of  physical  systems.  Physical  systems  are  almost 
always  described  by  a  system  of  differential  equations  and  in 
studvinq  the  behavior  of  a  physical  system,  we  ultimately  study  the 
behavior  of  the  solutions  to  the  differential  equations  which 
describe  the  system.  Thus  the  study  of  stability  of  a  nhysical 
system  is  reduced  to  a  consideration  of  the  stability  of  solutions 
of  differential  equations,  and  stability  theory  is  today  an  intearal 
part  of  that  segment  of  analysis  known  as  the  theory  of  ordinary 
differential  equations. 


Stability  theory  is  quite  advanced  today  and  there  are  numerous 
concents  and  approaches  which  vary  with  both  the  form  and  order 
of  the  system  of  differential  equations.  For  example,  the  Poincare- 
Bendixson  theory  is  an  analytical-toooloqical  approach  which  is 
Darticularly  useful  for  second  order  autonomous  systems,  while  for 
non-autonomous  systems,  the  concept  of  stability  in  the  sense  of 
Lyanunov  is  utilized  and  the  approach  is  purely  analytical.  This 
paper  is  concerned  with  a  general  study  of  stability  in  the  sense  of 
Lyaounov,  with  particular  emphasis  on  systems  of  linear  differential 
equations  though  some  results  for  nonlinear  systems  are  considered. 
Section  II  deals  with  a  review  of  familiar  facts  from  the  theory  of 
linear  ordinary  differential  equations.  Stability  and  boundedness  for 
linear  systems  are  discussed  in  section  III.  Non-homoqeneous  systems 
are  considered  in  section  IV,  and  the  concents  of  uniform  and  restric- 
tive stability  are  treated  in  section  V.  Some  results  for  non-linear 
systems  are  given  in  section  VI,  and  concluding  remarks  follow  in 
section  VII. 


II.   PRELIMINARY  REMARKS 


We  consider  the  behavior  of  solutions  of  the  nth  order  linear 

non-homoqeneous  system 

(2.1)  x1  =  A(t)x 

where  x  =  x(t)  is  an  n  by  1  column  vector  and  A(t)  is  an  n  by  n  matrix 

with  elements  a..(t),  i ,  j  =  1 , . . .  r\,   which  are  real  valued  functions 

continuous  on  the  half-line  [t  ,°°) . 

o 

The  norm  of  x,  A  is  defined  by 
(2.21)  || 

(2.211)  l|A(t)||  =  E|Au(t)| 

and  note  that  both  |x||  and  |A| |  are  functions  of  t.  The  usual 
oroDerties  of  norms  are  satisfied  and  are  stated  here  without  proof. 
For  |x|  they  are , 

I  x  |  |  L  0   and 
llxll  -  Myll 


x(t)||  =  z|x.(t) 


(2.31) 

(2.311) 

(2.3iii) 

(2.3iv) 

and  for  |A| 

(2.41) 

(2.411) 

(2.4iii) 

(2.4iv) 

(2.4v) 


aX| 
/x(t)dt 


<  J 


|x| |  =  0  iff  x(t)  =  0, 
<  |  |x±y|  |  <J  |x|  |  +  |  |y|  |  , 
x| |  for  a  a  real  number, 
x(t) I |dt, 


0   and  | |A| |  =  0   iff  A(t)  =   0, 


A] 


B 


A±B 


A 


AB| 
Ax 

aA 


A 


for  a  a  real  number. 


Canital  letters  will  denote  n  by  n  matrices,  lower  case  letters 
n  by  1  column  vectors,  and  Greek  letters  real  numbers  unless  otherwise 


noted.  Primes  (')  denote  differentiation  with  resnect  to  t;  A  denotes 
the  transpose  of  A;  tr(A)  denotes  the  trace  of  A;  det  A  denotes  the 
determinant  of  A;  and  adj  A  decotes  the  adjoint  of  A.  Lower  subscripts 
denote  components  of  a  vector  or  matrix,  for  example,  x.  is  the  ith 
component  of  the  vector  x. 

The  questions  of  existence  and  uniqueness  of  solutions  to  the 
system  (2.1)  may  be  dealt  with  quite  briefly.  Solutions  to  (2.1) 
satisfying  the  initial  condition  x(t  )  =  x  exist  and  are  unique  on 
any  interval  in  which  A(t)  is  continuous.  Hereafter,  continuity  of 
A(t)  on  [t  ,°°)  will  be  assumed  so  that  existence  and  uniqueness  of 
solutions  for  systems  of  this  form  are  always  assured. 

Linear  homogeneous  systems  of  differential  equations  are  charac- 
terized by  the  fact  that  any  linear  combination  of  solutions  is  also 
a  solution.  The  space  of  all  possible  solutions  to  the  system  (2.1) 
constitutes  an  n-dimensional  vector  space.  Thus  there  exists  n 
linearly  independent  solutions,  <j>  ,<)>,...,  <j>  ,  which  form  a  basis  for 
this  vector  space  of  solutions.  Any  particular  solution  x(t)  of  (2.1) 
can  then  be  expressed  as  a  linear  combination  of  elements  in  this 
basis.  Since  any  n  linearly  independent  solutions  constitute  a  basis, 
it  is  convenient  to  pick  a  particular  set  and  define  the  fundamental 
matrix  X(t)  as  follows. 

(2.5)  X(t)  =  4>j(t)   (i,j  =  l,...,n), 

i.e.,  the  columns  of  X(t)  are  n  linearly  independent  solutions  of 
(2.1)  and  these  n  solutions  ^    (i  =  l,...,n)  are  chosen  to  satisfy  the 

initial  condition  a1. (t  )  =  &..        (i.j  =  l,...,n)  so  that  X(t  )  =  I, 
Tj  o     ij  o 

the  identity  matrix. 


From  (2.5),  we  note  that  X(t)  satisfies  the  matrix  differential 
equation 

(2.6)  X'  =  A(t)X, 

and  that  any  solution  x(t;t  ,x  )  of  (2,1)  satisfying  the  initial 

condition  x(t  )  =  x  may  be  written 

o    o  J 

(2.7)  x(t)  =  X(t)xQ. 

From  the  formula  for  the  derivative  of  the  determinant  of  a 

matrix, 

_d_detX(t)  =  tr  [A(t)]  det  X(t), 
dt 

we  obtain  the  Jacobi-Liouville  formula 

t 

(2.8)  det  X(t)  =  det  X(t  )  exD  /  tr  A(s)  ds 

t 

=  exD  /  tr  A(s)  ds, 

t 
o 

and  note  from  (2.8)  that  X(t)  is  non-sinqular  for  t  >_  t  . 

Any  non-singular  matrix  W(t)  which  is  a  solution  of  (2.6)  mav 

serve  as  a  basis  matrix  though  if  W(t  )  f   I,  we  no  longer  have  (2.7) 


but  instead 


where  w  f   x 
o    o 


x(t;tQ,xo)  =  W(t)  wQ 


From  the  formula  for  the  derivative  of  the  inverse  of  a  matrix, 

d  [X_1(t)]  =  -X-](t)  X'(t)  X'](t), 
HI  ] 

we  obtain  the  result  that  if  Y(t)  =  X  (t)  ,  then  Y(t)  satisfies 

the  matrix  differential  eguation 

(2.9)  V     =  -A(t)T  Y, 


and  since  Y(t)  is  non-singular  and  Y(t  )  =  I,  we  note  that  Y(t) 
is  the  fundamental  matrix  to  the  adjoint  system  of  (2.1)  qiven  by 

(2.10)  x'  =  -A(t)T  x. 

The  fundamental  matrix  X(t)  Drovides  an  explicit  form  for  the 
solution  of  a  non-homogeneous  system  of  the  form 

(2.11)  x'  =  A(t)  x  +  b(t) 
as  the  following  proDosition  illustrates. 

Proposition  2.1 .  The  solution  x(t;t  ,x  )  of  the  system  (2.11), 
where  A(t)  and  b(t)  are  continuous  on  [t  >°°)  may  be  written 

x(t)  =  X(t)  x  +  /  X(t)  X_1(s)  b(s)ds, 

where  X(t)  is  the  fundamental  matrix  of  the  homogeneous  system 
x'  =  A(t)  x. 

Proof:  Direct  substitution. 

A  similar  formula  may  be  obtained  for  non-linear  mixed  systems 
of  the  form 

(2.12)  x'  =  A(t)  x  +  f(t,x). 

Proposition  2.2.  Let  f(t,x)  be  continuous  for  t  >  t  and 
—  o 

|x|  <  £,  and  A(t)  continuous  for  t  >_  t  .  Then  the  solution 
x(t;  t  ,x  )  of  (2.12)  may  be  written 

x(t)  =  X(t)  x  +  /  X(t)  X_1(s)  f(s,x(s))  ds, 

t0 

where  X(t)  is  the  fundamental  matrix  of  x'  =  A(t)  x. 

Proof:  Direct  substitution. 

Both  Propositions  2.1  and  2.2  simplify  in  the  case  A(t)  =  A, 
a  constant  matrix  and  t  =  0.  For  then,  the  matrix  X(t)  X"  (s) 
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satisfies  (2.6)  with  initial  condition 

X(t)  X"T(s)  =  I 
at  t  =  s,  and  the  same  follows  for  X(t-s).  Hence  by  uniqueness, 
(2.13)  X(t)  X_1(s)  =  X(t-s), 

for  A  a  constant  matrix  and  X(0)  =  I. 

Other  well  known  results  from  the  theory  of  ordinary  differential 
equations  will  be  stated  when  needed. 
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III.   LINEAR  HOMOGENEOUS  SYSTEMS 
Consider  a  general  first  order  system  of  ordinary  differential 
equations  of  the  form 
(3.1)  x'  =  f(t,x) 

where  x(t)  is  an  n  by  1  column  vector  and  f(t,x)  is  an  n  by  1 
column  vector  defined  and  continuous  in  the  region  I  =  {(t,x): 
0<_t<°°,   |x|  <  °°}. 

We  wish  to  define  a  concept  of  stability  for  this  system.  One 
definition  which  immediately  comes  to  mind  might  be  to  define  the 
system  to  be  stable  if  it  exhibits  only  finite  motions,  i.e.,  all 
solutions  remain  bounded  as  t->°°.  However,  such  a  conceDt  is  not 
particularly  useful  as  it  distinguishes  variational  effects  only  as 
to  their  being  finite  of  not.  For  example,  an  earthquake  would  be 
stable  in  this  sense. 

What  we  wish  to  embody  in  a  stability  concent  is  a  means  by  which 
we  can  measure  the  effects  of  small  variations  on  the  system  (or  a 
particular  solution  of  the  system).  One  criterion  we  would  desire  of 
a  stable  system  is  that  small  variations  not  affect  the  system  in  a 
drastic  manner.  This  is  the  nrimary  consideration  embodied  in  the 
concent  of  Lyaounov  stability,  which  imposes  stringent  restrictions 
on  the  solutions  of  (3.1),  by  requiring  that  solutions  which  are 
initially  close  together  remain  close  together  for  all  future  time. 
Formally,  we  define  stability  as  follows. 
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Definition  3.1 .  Let  x(t;t  ,x  )  be  a  solution  of  (3.1)  satisfyinq 
(3.2i)  x(t)  exists  for  t  >_  tQ, 

and 

(3.2ii)  (t,x(t))e  r   for  t  >  tQ. 

Then  x(t)  is  said  to  be  stable  (in  the  sense  of  Lyanunov)  if 
(3.2iii)         there  exists  a>0  such  that  every   solution 

y(t;tQ,yo)  exists  on  [t0,°°)  and  [t,y(t)) 


whenever  |x  -  y  I  <  a 

M  0   -  0 1  I 


and 


(3.2iv)  given  e>0  there  exists  6  =  6(e,f,x  ), 

0<6<a,  such  that  |  |x(t)  -  y(t)  |  <  £  for  t  >_  t 

for  all  solutions  y(t)  with  | | x  -  y  |  <  6. 

Strictly  speaking,  the  above  defines  stability  at  the  right,  since 

the  behavior  as  t+°°  is  considered;  replacinq  the  interval  [t  ,°°)  by 

the  interval  (-°°,t  ]  would  define  stability  at  the  left.  Note 

from  (3.2iv)  that  stability  is  equivalent  to  uniform  continuity  of 

x(t;t  ,x  )  on  [t  ,°°)  with  respect  to  the  initial  vector  x  .  A  solution 
ooo  o 

which  is  not  stable  is  said  to  be  unstable.  Hereafter,  stability  will 

always  refer  to  Lyapunov  stability  at  the  right  as  defined  above. 

A  special  case  of  stability  occurs  when  solutions  which  are 

initially  close  approach  one  another  as  t+«>.  This  leads  to  the 

followinq  definition. 

Definition  3.2.  A  solution  x(t;t  ,x  )  of  (3.1)  is  said  to  be 
v   o  o 

asymntotically  stable  on  [t  ,°°)  if  it  is  stable  and  in  addition  there 
exists  6'  >  0  such  that  |  |x(t)  -  y(t)||+0  as  t-*»  whenever  y(t;t  ,y  ) 
is  a  solution  with  llx  -  y  I  <  6'. 

ll  0   JQ\\ 
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Note  that  by  definition,  a  solution  which  is  asymptotically  stable 
is  stable.  If  all  solutions  of  the  system  (3.1)  are  (asymntotically) 
stable,  then  we  say  that  the  system  is  (asymptotically)  stable,  while 
if  there  exists  an  unstable  solution  to  (3.1),  we  will  call  the  system 
unstable. 

The  following  examples  illustrate  some  of  these  concepts. 

Examnle  3.1  .  Consider  the  one-dimensional  system 

on  [0,oo).  The  general  solution  x(t;0,x  )  is  given  by  x(t)  =  x  e"  , 

so  that 

I |x(t)  -  y(t)| I  =  |x  -  v  I  e_t  <  |x  -  y  I 
I  I  v  '   J\'i.i     '  o   -  o1     —  '  0    o1 

Conditions  (3.2  i,ii,  and  iii)  are  satisfied  as  is  (3.2iv)  with 

6  =  e.  Thus  the  system  is  stable. 

In  addition, 

I  |x(t)  -  y(t)|  I  =  |x  -  y  I  e_t  ■*  0  as  t-» 

so  that  the  system  is  asymptotically  stable. 

Example  3.2.  Every  non-trivial  solution  of  the  one-dimensional 

system 

x'  =  0   on  [0,oo) 

is  of  the  form  x(t;0,x  )  =  x  and  thus 

0    0 

I |x(t)  -  y(t)| I  =  |xQ  -  yQ\  . 
Conditions  (3.2i,ii,  and  iii)  are  satisfied  as  is  (3.2iv),  aaain  with 
6  =  e. 

However, 

lim  I  I x ( t)  -  y(t) I   =  |x  -  v  I  >  0, 

so  that  no  solution  is  asymptotically  stable.  Thus  the  system  is 
stable. 
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Example  3.3.  Non-trivial  solutions  of  the  one-dimensional  system 

x1  =  x   on  [0,°°) 

are  of  the  form  x(t;0,x  )  =  x  e  and  thus 

o      o 

I  |x(t)  -  y(t)|  I  =  |x"  -  y  le* 
which  can  be  made  arbitrarily  larqe  independent  of  |x  -  y  |  by  takinq 
t  large  enough.  Thus  the  system  is  unstable. 

Example  3.4.  Cesari  [2]  shows  that  the  trivial  solution  z(t)  =  0 


of  the  one  dimensional  system 

.2 


on  [0,°°) 


is  not  stable  since  every  (real)   solution  x(t;0,x   )   is  of  the  form 

x(t)     =        xo 

XPExp 

and  solutions  with  x  >  0  cease  to  exist  at  t  =  (x  )"  in  violation 
of  (3.2iii). 

For  linear  systems  of  the  form 
(3.3)  x'  =  A(t)  x  +  b(t) 

with  A(t)  and  b(t)  continuous  on  [t  ,°°),  (3.2i,ii,  and  iii)  are  always 
satisfied  so  that  (3.2iv)  becomes  the  determining  factor  as  to  whether 
or  not  solutions  to  these  systems  are  stable.  For  these  systems, 
we  have  the  fol^jwing  useful  proposition. 

Proposition  3.1 .  If  the  linear  system  (3.3)  has  a  stable  solution 

then  all  its  solutions  are  stable,  i.e.,  the  system  is  stable. 

Proof:  Let  x(t;t  ,x  )  be  a  stable  solution  and  y(t;t  ,y  ) 
00  ^ v   o  o 

be  any  other  solution.  Let  e  >  0. 

If  z(t;t  ,z  )  is  a  solution  of  (3.3),  then  x(t)  -  y(t)  +  z(t) 
is  also  a  solution  since  the  system  is  linear.  Then  since  x(t) 
is  stable,  there  exists  6  >  0  such  that  if 
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||x(t  )  -  [x(t  )  -  y(t  )  +  z(t  )]|   =   ly  -  z 

11    0        0    J v  0      v  0   ' '      M Jo    0 

<  6 

then 

||x(t)  -  [x(t)  -  y(t)  +  z(t)]||  =  ||y(t)  -  z(t) 

<  e 
which  says  that  y(t)  is  stable.  Thus  all  solutions  are  stable. 

Corollary  1 .  Proposition  3.1  remains  true  if  stable  is  reDlaced 
by  asymptotically  stable. 

Proof:  The  proof  is  identical  to  that  above. 

We  next  define  the  concept  of  boundedness  of  solutions  to  3.1. 

Definition  3.3.  A  solution  x(t)  of  (3.1)  is  said  to  be  bounded 
if  there  exists  M  >  0  such  that    x ( t ) I  <  M,  t  >  t  , 
i.e.,  if  ||x(t)|  is  a  bounded  function  of  t  on  [t  ,'»).  If  all 
solutions  to  (3.1)  are  bounded,  we  say  that  the  system  is  bounded. 

Stability  and  boundedness  are  independent  concepts  as  the 
following  examDles  illustrate. 

Example  3.5.  All  solutions  of  the  one  dimensional  system 

x1  =  2t   on  [O.oo) 

2 

are  of  the  form  x(t;t  ,x  )  =  x  +  t  and  are  stable  yet  unbounded 

oo      o 

on  [0,°°). 

Example  3.6.  Cesari  [2]  shows  that  every   non-trivial  solution  of 
the  two  dimensional  system 

x'  =  y 

y'  =  -  1/2  x  (x2  +  (x4  +  4y2)1/2s 
or  equivalently, 

z"  =  -1/2  z  (z2  +  (z4  +  4(z')2)1/2 
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is  of  the  form 

z(t)  =  c  sin  (ct+d),   c  f     0 
and  is  thus  bounded;  none  of  these  solutions  is  stable. 

For  linear  homoqeneous  systems  of  the  form 
(3.4)  x'  =  A(t)  x 

stability  and  boundedness  are  equivalent  as  the  followinq  theorem 
illustrates. 

Theorem  3.1 .  The  system  (3.4)  is  stable  iff  it  is  bounded. 

Proof:  Assume  the  system  is  stable.  To  show  that  it  is  bounded, 

it  suffices  to  show  that  there  exists  M  >  0  such  that   |X(t)|  <^  M 

where  X(t)  is  the  fundamental  matrix  of  the  system.  Since  the  system 

is  stable,  in  particular,  the  trivial  solution  z(t)  =  0  is  stable. 

Thus  if  e  >  0  there  exists  &  >  0   such  that  every  solution  x(t;t  ,x  ) 

j  v   o  o 

with  ||x  ||  <  6  satisfies  ||x(t)||  <  e   for  t^t  .  In  particular,  let 

x(t)  be  that  solution  with  initial  vector  x  given  bv 

o 

x   =  f  e1 

o    2 

where  e  is  the  vector  whose  jth  component  is  qiven  by  e.  =  6... 

Then  |x    <  6  so  that 

ii  0i  i 

||x(t)||   =   ||X(t)  x0||  <  e,   t  >  tQ. 

r\  i  "i 

But  for  this  choice  of  x  ,  X(t)  x   =  t  <f>  ,  where  d>  is  the  ith  column 

oo    c 


of  X(t).  Thus 


6_  |  I*""  (t)  |  I  <  e   t  >  t 

2  "  ° 


This  procedure  is  valid  for  i  =  l,...,n  so  that 

l|X(t)||  -   I  ||*1(t)|| 


1=1  6 


and  hence  X(t)  is  bounded.  Thus  the  system  is  bounded 


17 


Conversely,  assume  the  system  is  bounded  so  that  X(t)  <^  M, 
for  t  >_  t  .  Then  if  x(t;t  ,x  )  and  y(t;t  ,y  )  are  two  solutions  to 
(3.4) 

||x(t)  -  y(t)||  =  ||X(t)  (xq  -  yo)||  <  M  ||xq  -  yQ 
so  that  if  |  |xQ  -  yQ|  |  <  jjj-  ,  then  |  |x(t)  -  y(t)|  |  <  e   for  t  >_  t 
which  shows  that  the  system  is  stable. 

Corollary  1 .  The  system  x1  =  A(t)  x  is  asymototically  stable 
iff  X(t)->0  as  t+°°,  where  X(t)  is  the  fundamental  matrix  of  the  system. 

Proof:  The  reverse  implication  is  obvious.  To  show  the  forward 
implication,  note  that  the  hypotheses  imply  that  the  trivial  solution 
is  asymptotically  stable  and  the  result  then  follows  as  in  the  proof 
of  the  first  half  of  Theorem  3.1. 

Before  proceeding  to  consider  various  results  concerning  stability 
in  linear  homogeneous  systems,  it  is  convenient  to  introduce  the 
following  lemma,  known  as  Gronwall's  Lemma,  which  is  an  extremely 
useful  tool  in  stability  theory. 

Lemma  3.1.  If  f(t)  and  g(t)  are  real  valued  non-neqative  functions, 

f(t)  is  continuous  on  [t  ,°°)  and  g(t)  is  integrable  in  every   finite 

interval  contained  in  [t  ,»)  and  if  for  some  constant  a  >  0 

o       .  — 

(3.5)  f(t)  <  a  +     f     f(s)  g(s)  ds, 

t 
o 

then 

t 

(3.6)  f(t)  ■:  a  exp  /   g(s)  ds. 

t 
o 

Proof:  If  a  >  0,  then  we  have 

f(t)  g(t)      ig(t). 

a  +  /*  f(s)g(s)ds 
t 
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Integrating  the  above  between  the  limits  t  and  t  yields 

t  t 

In  [a  +  /  f(s)g(s)ds]  -  In  a  <  /  g(s)ds, 

t  t 

o  o 

or  taking  exDonents, 

t  t 

a  +  /  f(s)g(s)ds  <  a  exp   /  g(s)ds, 

t  t 

o  o 

and  thus 

t 
f (t)  <  a  exp   /  g(s)ds. 

If  a  =  0,  then  (3.5)  holds  for  all  real  numbers  B  >  0,  and  thus 
(3.6)  holds  which  implies  as  b+0  that  f(t)  =  0.  This  completes  the 
Droof . 

Utilizing  Lemma  3.1,  it  is  Dossible  to  Drove  the  following 
proDosition. 

00 

Proposition  3.2.  If   /||A(t)|  dt  <  °°,  then  the  system  x'  =  A(t)  x 

is  stable. 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  system.  Then 

from  (2.6) 

t 
X(t)  =  I  +  /  A(s)  X(s)  ds, 

so  that 

t 

||X(t)||  :  n  +  /  ||A(s)||  ||X(s)||  ds. 
t 

0 

Applvinq  Gronwall's  Lemma  yields 

t 

|  |X(t)|  |  <  n  exD  /  |  |A(s)  |  |  ds  <_  n  exp  /  |  |A(s)|  |  ds 

t  t 

o  o 

and  thus  X(t)  is  bounded.  By  Theorem  3.1,  the  system  is  stable. 
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We  next  consider  the  linear  homogeneous  system  with  constant 
coefficients 

(3.7)  x'  =  A  x 

whose  fundamental  matrix  X(t)  is  given  by 

(3.8)  X(t)  =  exp  [(t-tQ)A], 
where  exn  (A)  is  defined  by 

exp  (A)  =  I  +   i  1,  An  . 
n=ln* 

The  stability  of  the  system  (3.7)  is  directly  related  to  the  nature 

of  the  roots  of  the  characteristic  polynomial  of  A  given  by 

(3.9)  det  (A  -  Al)  =  0. 

Since  (3.9)  is  an  nth  degree  polynomial  equation  in  a,  there  exists 
n  roots  to  this  equation  (including  repeated  roots)  which  may  be 
real  or  complex.  Corresponding  to  each  root  is  a  solution  of  (3.7) 
and  the  n  roots  thus  furnish  n  linearly  independent  solutions  to 
(3.7).  If  a  is  a  real  root  of  multiplicity  k,  then  corresponding 

to  a  are  the  k  linearly  independent  solutions 

nxt     .At      .k-1  At 
e  ;  te  ;  . . . ;  t   e 

while  if  a+i'm  is  a  complex  root  of  multiplicity  m,  then  so  is 
A-iy  and  corresponding  to  these  2m  roots  are  the  2m  linearly  indepen- 
dent solutions 

„At    .   At  .   .   .  At    .   .  At  .   ,      jn-1     .   .m-1  .   . 
e  cos  yt,  e  sin  yt;  te  cos  yt,  te  sin  yt;  . . . ;  t   cos  yt,  t   sin  yt, 

Since  sin  t  and  cos  t  are  bounded,  and  since  for  a,  6  >0, 

lim  tae    =  0,  we  note  that  if  the  real  parts  of  the  eigenvalues  of 

t-><» 

A  (roots  to  (3.9))  are  negative,  then  all  solutions  of  (3.7)  will 
converge  to  zero.  When  this  is  the  case,  we  have  the  following 
definition. 
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Definition  3.4.  If  the  real  parts  of  the  eiqenvalues  of  A  are 
negative,  then  we  say  that  the  characteristic  nolynomial  of  A  is 
stable. 

When  the  characteristic  polynomial  of  A  is  stable,  then  the 
solutions  of  (3,7)  are  of  exponential  order  as  the  following  lemma 
illustrates. 

Lemma  3.2.  If  the  characteristic  nolynomial  of  A  if  stable, 
then  all  solutions  of  (3.7)  satisfy 
(3.10)  ||x(t)||  1Me"at   t  >  tQ, 

where  M  and  a  are  positive  constants. 

Proof:  Let  A,   =  y.  +  io,  ,   k  =  1 ,. . .  ,m,  m  <_  n,  be  the 

distinct  eigenvalues  of  A.  Since  the  characteristic  nolynomial  of  A 
is  stable,  there  exists  a  >  0  such  that  y,  +  a  <  0,  k  =  l,...m. 
A  solution  z. (t)  corresponding  to  x.    is  of  the  form 

zk(t)  =  tr  eAkt  bk, 
where  b,  is  a  constant  vector  having  the  entry  1  in  one  of  its  rows 
and  zeros  elsewhere,  and  where  r  depends  on  the  multiplicity  of  A.  . 


Then  since  p,  +  a  <  0, 


||zk(t)||  eat>0   as  t- 


and  thus  there  exist  constants  M,  such  that 

If  x(t)  is  any  solution  to  (3.7),  it  can  be  expressed  .as  the  linear 
combination 

x(t)  =   EBi  z^t), 
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where  the  summation  is  over  the  n  basis  solutions.  Then 

||x(t)||  1E|3.|  ||z.(t)|| 

<  max  1 6. 1  EM.  e'at  <  M  e~a1\ 
—     I  M-|  i     t      _       » 

where  M  =  max  |  B-  |  £M- . 

It  is  now  easy  to  show  that  when  the  characteristic  Dolynomial 
of  A  is  stable,  then  (3.7)  is  an  asymptotically  stable  system. 

Proposition  3.3.  If  the  characteristic  polynomial  of  A  is  stable, 
then  the  system  x'  =  Ax  is  asymptotically  stable. 

Proof:  By  Lemma  3.2,  every  solution  x(t)  of  (3.7)  satisfies 
|x(t)|  <_  M  e~a   for  appropriate  M  and  a.  Since  the  columns  of 
X(t)  are  solutions,  where  X(t)  is  the  fundamental  matrix  of  (3.7), 
we  have 

I |X(t)| I  <  Mne"at   t  >  t  , 
and  thus  |X(t)||-»-0   as  t-*»  and  the  system  is  asymptotically  stable 
by  Corollary  1  to  Theorem  3.1. 

Thus  the  stability  of  a  linear  homogeneous  system  with  constant 
coefficients  depends  solely  upon  the  nature  of  the  eigenvalues  of  A. 
The  question  now  arises  as  to  the  stability  of  a  system  (3.4)  in  which 
A(t)  is  not  constant.  Proposition  3.1  insures  that  inteqrability  of 
A(t)  is  sufficient  to  guarantee  stability  in  the  system.  However,  this 
is  a  rather  stringent  reguirement  and  is  of  little  practical  signi- 
ficance in  determining  the  stability  of  a  given  system.  It  is  often 
easier  to  determine  stability  in  a  given  system  by  comparing  it  to  a 
known  stable  system.  The  simplest  case  in  which  this  may  be  done 
successfully  is  when  A(t)  is  asymptotic  to  a  constant  matrix  B. 
The  following  theorem  then  applies 
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Theorem  3.2.  If  the  characteristic  polynomial  of  B  is  stable, 
A(t)  is  continuous  on  [0,°°),  and 

oo 

/  |  |A(t)  -  B|  |  dt  <  °°  , 
0 

then  the  system  x'  =  A(t)  x  is  asymptotically  stable. 

Proof:  The  system  x'  =  A(t)  x  may  be  written  in  the  form 

xf  =  Bx  +  [A(t)  -  B]  x 

and  in  this  manner  may  be  treated  as  though  it  were  a  non-homogeneous 

system  with  constant  coefficients.  The  solution  x(t;0,x  )  may  then 

be  written  by  Proposition  2.2  and  (2.13) 

t 
x(t)  =  X(t)  x  +  /  X(t-s)  [A(s)  -  B]  x(s)  ds, 
0    0 

where  X(t)  is  the  fundamental  matrix  of  the  constant  system 

x1  =  B  x.  By  Lemma  3.2.,  there  exist  M,a  >0  such  that 

|  |X(t)|  |  <_  Me"at   and  thus 

||x(t)||  eat1M  ||x  ||  +  m   ||A(s)-B||  ||x(s)||  eaSds 

o 

Now  applying  Gronwalls'  Lemma, 

t  t 

|  |x(t)|  |  eal1M  |  |x  |  |  exp  /  M  |  |A(s)-B|  |  ds, 

o 
from  which 

|x(t)|  eaL  <  M  |x  |  exn  /  M  |A(s)-B| ;  ds 

0      0 

i*. 

where  K  is  a  convenient  constant.  Hence  |x(t)|  <  Ke"a   and  if 
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y(t;0,y  )  is  any  other  solution, 

I | x(t)  -  y(t)| |  <  2Ke"at  +  0   as  t*» 
and  thus  the  system  x'  =  A(t)  x  is  asymptotically  stable. 

Theorem  3.2.  may  be  modified  in  the  event  A(t)  does  not  differ 
much  from  B,  i.e.,   |A(t)  -  B| |   is  small  as  follows. 

Theorem  3.3.  If  the  characteristic  Dolynomial  of  B  is  stable, 
and  |A(t)  -  B| |  c  on  [0,»)  for  sufficiently  small  c,  then  the 
system  x1  =  A(t)  x  is  asymptotically  stable. 

Proof:  Let  x(t;0,x  )  be  a  solution  of  x'  =  A(t)  x.  Then  x(t) 

may  be  written 

•   t 
x(t)  =  X(t)  xn  +  /  X(t-s)[A(s)  -  B]  x(s)  dss 
0   0 

where  X(t)  is  the  fundamental  matrix  of  x '  =  B  x.  Then  proceeding  as 

in  the  proof  of  Theorem  3.2,  we  have  as  before 

t  t 

|  |x(t)|  |  eaL  <_  M  |  |x  |  |  exp  M  /  |  |A(s)  -  B|  |  dss 

o 
where  M  and  a  are  as  previously  defined  in  Theorem  3.2.  Thus 

||x(t)||  <_  M  ||x  ||  exp  (Mc-a)  t  , 
so  that  if  c  is  sufficiently  small,  i.e.,  c  <  «•»  then  x(t)->0  as 
t->oo  and  since  x(t)  is  an  arbitrary  solution,  the  system  is 
asymptotically  stable. 

If  the  characteristic  polynomial  of  B  is  not  stable  but  the  solu- 
tions of  x'  =  B  x  are  bounded,  then  the  following  theorem  holds. 

Theorem  3.4.  If  all  solutions  of  x'  =  B  x  are  bounded  on  [0,«°), 
A(t)  is  continuous  on  [0,°°)  and 

oo 

/  I  |A(t)  -  B|  |  dt  <  °°s 
0 

then  the  system  x'  =  A(t)  x  is  stable. 
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Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  constant 
system  x'  =  B  x.  Then  there  exists  M  >  0  such  that  |X(t)|  <_  M 

for  t  >_  0.  If  x(t;0,x  )  is  a  solution  of  x'  =  A(t)  x,  then 

t 
x(t)  =  X(t)  x  +  /  X(t-s)[A(s)  -  B]  x(s)  ds. 
0    0 


Thus 


t 

l|x(t).||  cM  ||xo||  +   /M  ||A(s)  -  B||  ||x(s)||  ds, 


0 


and  applying  Gronwall's  Lemma 


t 
x(t)|  L  :  M  |  [x  J  |  exp  /  M  |  |A(s)  -  B|  |  ds, 
0       0 


from  which 


|x(t)|  |  :  M  |  |x  |  exp  /  M  |A(s)  -  B|  ||  ds 
0       0 

£  K  <  «> 

where  K  is  a  suitable  constant.  The  above  shows  that  all  solutions  are 
bounded  and  hence  the  system  x1  =  A(t)  x  is  stable  by  Theorem  3.1. 
Theorems  3.2,  3.3,  and  3.4  insure  that  if  |A(t)  -  B|  is  well 
behaved,  then  the  solutions  of  x'  =  A(t)  x  will  resemble  those  of  the 
constant  system  x'  =  B  x  as  t^°°  when  A(t)  is  asymptotic  to  B.  In 
view  of  this,  we  might  expect  these  results  to  generalize  to  the  case 
when  B  is  not  constant.  This  is  not  always  the  case  though  it  is 
possible  to  effect  a  comparison  between  two  non-constant  systems  under 
more  restrictive  hypotheses.  Toward  this  end,  we  Drove  the  following 
1 emma . 
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Lemma  3.3.  If  the  system  x*  =  A(t)  x  is  stable  on  [t  ,~) 

and 

t 

"Mm  inf  /  trA(s)  ds  >  -  °°  , 

t->°°      t 
o 

then  X"  (t)  is  bounded  on  [t  ,*»),  where  X(t)  is  the  fundamental  matrix 
of  the  system. 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  system. 
Then 

v-ht)     -  adJ  X(t) 
X  {t}     ~     det  X(t) 

and  since  x1  =  A(t)  x  is  stable,  it  follows  that  X(t)  is  bounded  for 

t  >  t  and  this  implies  that  there  exists  L  >  0  such  that 
—  o         K 

I  I ad j  X(t)| I  <  L   for  t  >  t  . 

1  '       '  '  —        —  o 

Then  this  relation  and  (2.8)  yield 

I'lx^tJll  =  — 


t 

exD  /  tr  A(s)  ds 
t. 


and  since 


t 
lim  inf  /  tr  A(s)  ds  >  -  °°, 
t+°°     t 


it  follows  that 


so  that 


exp   /  tr  A(s)  ds    k  >  0, 


X""1  Ct>  I  |  <..{=■  for  t  >_t0 


and  thus  X"  (t)  is  bounded. 
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We  can  now  state  the  following  comoarison  theorem  for  two  non- 
constant  systems. 

Theorem  3.5.  If  the  system  x'  =  B(t)  x  is  stable  on  [t  ,°°), 

t 

lim  inf  /  tr  B(s)  ds  >  -  °°, 

t-~»     t 
o 

A(t)  is  continuous  on  [t  ,°°)  and 

oo 

I      ||A(t)  -  B(t)||  dt  <  -, 

then  the  system  x'  =  A(t)  x  is  stable. 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  system 
x'  =  B(t)  x.  Then  the  hypotheses  imply  that  there  exists  M  >  0 
such  that  for  t"£t;  |  |X(t)|  |  <_M,  and  |  |X_1(t)|  |  <M. 
Then  if  x(t;t  ,x  )  is  any  solution  to  x'  =  A(t)  x,  x(t)  may  be  written 

1  1 

x(t)  =  X(t)  x  +  /  X(t)  X  '(s)  [A(s)  -  B(s)]x(s)ds, 

and  utilizing  the  above  we  obtain 

t  ? 
||x(t)||  <_M  ||x_||  +  /  V\l   ||A(s)  -  B(s)||  ||xCs]||ds 

Applying  Gronwall's  Lemma  yields 

t  - 
l|x(t)||  iM  ||x  ||  exp  /  M^  ||A(s)  -  B(s)|  ds , 


o 


lo 


so  that 


x(t)| |  <  M  | |x  I |  exD  /  M2  | |A(s)  -  B(s)| Ids 


o 


% 


where  K  is  a  convenient  constant.  Thus  all  solutions  of  x'  =  A(t)  x 
are  bounded  and  hence  stable  by  Theorem  3.1. 
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Theorem  3.5  gives  one  set  of  conditions  under  which  a  comparison 
may  be  effected  between  two  non-constant  linear  homoqeneous  systems 
and  is  a  well  known  result.  Another  set  of  conditions  which  have 
been  formulated  by  the  author  and  which  are  also  sufficient  to  effect 
a  comparison  are  stated  below  in  the  following  theorem. 

Theorem  3.6.  If  the  system  x'  =  B(t)  x  is  stable  on  [t  ,«>), 
A(t)  is  continuous  on  [t  ,°°)  and 

/  ||X-1  (t)  [A(t)  -  B(t)]  X(t)||  dt  <  - 

t 
o 

where  X(t)  is  the  fundamental  matrix  of  x'  =  B(t)  x,  then  the 
system  x'  =  A(t)  x  is  stable  on  [t  ,»). 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  x'  =  B(t)  x.  Then 
if  x(t;t  ,x  )  is  any  solution  to  x'  =  A(t)  x,  x(t)  may  be  written 

*       1 
x(t)  =  X(t)  xn  +  /  X(t)  X"'(s)  [A(s)-B(s)]x(s)  ds. 


From  this  we  obtain 


1  t   1 

X  '(t)x(t)  =  x  -.+  /  X"'(s)[A(s)  -  B(s)]x(s)  ds 

=  x  +  /  X-1(s)[A(s)  -  B(s)]X(s)X-1(s) 
0  x(s)ds 


and  thus 


X"1  (t)x(t)  |  |  <_  ||x0 


t 

t 


|X-1(s)[A(s)-B(s)]X(s)||  |  | X"1 (s) 


o 


x(s)  I  I  ds, 
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so  that  applying  Gronwall's  Lemma  yields 


-1  t     1 

X  '(t)x(t)||  <  ||x  ||  exp  /   ||X  '(s)[A(s)-B(s)] 

X(s)  I  ds. 


Therefore, 


X-1(t)x(t)||  <  ||x  ||  exp  /  | |X-1(s)[A(s)-B(s)] 


X(s)||  ds 


,-1 


1  K  < 
which  shows  that  X"'(t)x(t)  is  bounded. 

However, 


||x(t)||   =  ||  X(t)X-1(t)x(t)||  <  ||X(t)||  1 IX-1 (t) 

x(t)|| 
so  that  | |x(t)| |  <  K  | |X(t)| |  and  since  x'  =  B(t)  x  is  stable,  X(t)  is 
bounded  and  thus  x(t)  is  bounded.  Hence  all  solutions  to  x'  =  A(t)  x 
are  bounded  and  thus  stable  by  Theorem  3.1. 

Theorems  3.2  through  3.5  provide  conditions  under  which  the 
stability  of  a  given  system  may  be  determined  by  a  comnarison  with 
a  known  stable  system.  Other  conditions  under  which  similar 
theorems  can  be  established  have  been  given  by  Bellman  [1],  and 
Levinson  [6],  as  well  as  numerous  other  authors. 
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IV.  NON-HOMOGENEOUS  LINEAR  SYSTEMS 

The  systems  considered  here  are  of  the  form 
(4.1)  x'  =  A(t)  x  +  b(t) 

where  A(t)  and  b(t)  are  continuous  on  [t  ,°°). 

Recall  that  for  the  homogeneous  system,  b(t)  =  0,  boundedness 
and  stability  of  the  system  are  equivalent.  This  is  no  lonner  true 
for  the  system  (4.1)  as  the  following  example  illustrates. 

Example  4.1 .  The  one-dimensional  system 
x'  =  2t   on  CO,-) 
has  solutions  x(t;0,x  )  of  the  form 

x(t)  =  xQ  +  t2 
which  are  stable  yet  unbounded. 

The  above  example  shows  that  it  is  possible  for  a  non-homogeneous 
to  be  stable  but  unbounded.  The  reverse  behavior,  however,  is  not 
possible  as  the  following  theorem  shows. 

Theorem  4.1 .  If  the  system  x'  =  A(t)  x  +  b(t),  A(t)  and  b(t) 

continuous  on  [t  ,°°),  is  bounded  on  [t  ,°°)  then  it  is  stable  on 

1  0  0 

[t  ,-). 
L  0 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  homogeneous 
system  x'  =  A(t)  x.  If  (4.1)  is  bounded,  then  in  particular,  the 
solution  x(t;t  ,0)  given  by 

*   1 
x(t)  =  X(t)  /  X  '(s)  b(s)  ds 

*o 

is  bounded.  Thus  there  exists  Mn>  0  such  that  |x(t)|  <  M, ,  t  >  t  . 

11  v  ' ' '  —  r  —  0 
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If  y(t;t  ,y0)  is  any  other  solution  to  (4.1),  then 

t  1 

y(t)  =  X(t)yn  +  /  X(t)X"'(s)  b(s)  ds 

and  there  exists  M2  such  that 

||y(t)||  =  ||X(t)  yQ  +  x(t)||  <.M2  for  t  >_  tQ. 
Then 

||X(t)  yo||  -  ||x(t)||  iM2, 
so  that 

||X(t)  y0||  <H   t  LtQi 

where  M  =  M,  +  M0.  Therefore  if  w(t;t  ,w  )  and  v(t;t  ,v  )  are  any 
I    i  oo        oo 

two  solutions  of  (4.1)  and  if  e  >  0,  then  for  t  >_  t 


||w(t)  -  v(t)||  =  ||X(t)  (w0-vQ)  |  |  iM  ||w0-v0|| 


so  that  if  l|wa-v  ||  <  j|,  then  | |w(t)-v(t) | |  <e       t  >_  t, 
which  shows  that  (4.1)  is  stable. 

Theorem  4.1  admits  the  following  partial  converse. 

Theorem  4.2.  If  the  system  (4.1)  is  stable  and  one  solution  is 
bounded,  then  all  solutions  are  bounded. 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  homoqeneous 
system  x1  =  A(t)  x.  Once  again  the  first  part  of  the  nroof  consists 
of  showing  that  X(t)  is  bounded. 

If  (4.1)  is  stable,  then  in  particular,  the  solution  x(t;t  ,0) 
given  by 

1  1 

x(t)  =  /  X(t)X"'(s)b(s)  ds 

is  stable.  Therefore,  if  e  >  0,  there  exists  6  >  0  such  that 
for  t  >  t  , 
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! |w(t)-x(t)| I   =   I |X(t)w0| |  <  e 
whenever  w(t;t  ,w  )  is  a  solution  with  |w    <  6.  Lettinq 

V  g  0   0  M0M  3 

w  =  j  e1 ,  where  e  is  the  vector  whose  jth  comnonent  is  el  =  6..,  we 
obtain  from  the  above 


'Hi) 


<  2e_ 

6 

where  <j>  (t)  is  the  ith  column  of  X(t).  Then  since  this  is  valid  for 
i  =  1 ,...,n, 

||X(t)||  <  ^-    =   M. 


Now  let  z(t;t  ,z  )  be  a  bounded  solution  and  y(t;t  ,y  )  any  arbitrary 

0   0  0   0 

solution  of  (4.1).  Then  since  y(t)  =  z(t)  +  [y(t)  -  z(t)], 

l|y(t)||    ||z(t)||  +  ||X(t)[y  -z  ]|| 

i  i  -J    »   *  i  i   ,    ii    *   *  I  I       II    *   '  i"J  s*i  o   '  ' 

<| |z(t)|  I  +  M  I ly  -  z  I  I 
which  shows  that  y(t)  is  bounded.  Hence  the  system  (4.1)  is  bounded. 

Theorems  4.1  and  4.2  show  that  in  a  non-homogeneous  linear  system, 
boundedness  is  a  stronger  property  than  stability.  For  this  reason, 
it  is  usually  easier  to  determine  stability  in  a  system  (4.1)  than  it 
is  boundedness.  In  fact,  the  following  proposition  shows  that  the 
stability  of  (4.1)  may  be  determined  by  simDly  examing  its  homo- 
geneous part. 

Proposition  4.1 .  If  the  system  x'  =  A(t)  x  is  (asymptotically) 

stable  on  [t  »°°)»  then  so  is  the  system  x'  =  A(t)  x  +  b(t),  where 

b(t)  is  continuous  on  [t  ,°°). 

o 

Proof:  If  x(t;t  ,x  )  and  y(t;t  ,y  )  are  two  solutions  of 
oo         oo 

x'  =  A(t)  x  +  b(t),  then  the  difference  [x(t)-y(t)]  satisfies  the 
homogeneous  system  x1  =  A(t)  x.  Thus  since  the  trivial  solution 
z(t)  =  0  of  the  homogeneous  system  is  stable,  given  e  >   0,  there 


32 


exists  6  >  0  such  that  if  I (x  -  y  )  -  01 I  <  6  then 

1  '   0     0      '  ' 
||x(t)  -  y(t)  -  0||  <  e   t  £t0, 

which  is  precisely  the  statement  that  the  non-homogeneous  system  is 

stable.  The  proof  is  the  same  if  the  homogeneous  system  is,  in 

addition,  asymptotically  stable. 

If  A(t)  is  a  constant  matrix  A,  then  (4.1)  becomes 
(4.2)  x'  =  A  x  +  b(t) 

and  boundedness  is  easier  to  determine.  If  the  characteristic 
polynomial  of  A  is  stable,  then  the  following  theorem  holds. 

Theorem  4.3.  If  the  characteristic  polynomial  of  A  is  stable, 
b(t)  is  continuous  and  |b(t)|  <_  K,  for  t  >_  0,  then  the  system  (4.2) 
is  bounded. 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  constant 
system  x'  =  A  x.  By  Lemma  3,2,  there  exist  M  and  a  such  that 
|  |X(t)  |  |  <  Me"at  for  t  >_  0.  If  x(t;t  ,x  )  is  a  solution  of  (4.2) 


then  from 


t 
x(t)  =  X(t)x  +  /  X(t-s)b(s)  ds, 
0    0 


we  deduce  successively 

I  |x(t)| |    <  M   | |x 


o 


+       /  Me"a(t"s)K  d; 


1M   Hi   I'   +^<(l-e-at) 

O  a 

:M    i|.X    ||+M 

—        M    o1  '  a 

which  shows  that  x(t)  is  bounded.  Thus  the  system  (4.2)  is  bounded 

A  similar  theorem  holds  when  the  solutions  of  x'  =  A  x  are 
bounded,  though  the  forcing  term,  b(t),  must  be  integrable. 
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Theorem  4.4.  If  the  constant  system  x'  =  A  x  is  bounded  on 
[0,°°),  b(t)  is  continuous  on  [0,°°)  and 

00 

/  I |b(t)| |  dt  <  -  , 
0 

then  the  system  x'  =  A  x  +  b(t)  is  bounded. 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  x1  =  A  x.  Then 

there  exists  M  >  0  such  that  |  |X(t)  |  |  <_  M   t  >_  0.  If  x(t;0,x  )  is 

a  solution  to  x'  =  A  x  +  b(t),  then  from 

t 
x(t)  =  X(t)  x  +  /  X(t-s)b(s)  ds, 


0    0 


we  deduce  successively 


t 

I |x(t)| |  <  M  | |x  | |  +  /  M  | |b(s)| |  ds 
0     0 

oo 

_  M  I |x  I  I  +  /  M  | |b(s) | |  ds 
0     0 

<_  K  <  °° 

where  K  is  a  convenient  constant.  Thus  the  system  x'  =  A  x  +  b(t) 

is  bounded. 

The  boundedness  of  a  non-homoqeneous  system  x1  =  B(t)  x  +  b(t) 
when  B(t)  is  asymptotic  to  a  constant  matrix  A  may  be  determined  by 
effecting  a  comparison  between  it  and  the  constant  svstem  x'  =  A  x. 
The  following  theorem  gives  one  instance  when  this  may  be  done 
successfully. 

Theorem  4.5.  If  the  characteristic  polynomial  of  A  is  stable, 

b(t)  is  continuous  and  bounded  on  [0,°°),  B(t)  is  continuous  on  [0,°°) 

and  r   |B(t)  -  A|  dt  <  °°,  then  the  system  x'  =  B(t)  x  +  b(t)  is 

0 
bounded  and  asymDtotical ly  stable. 
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Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  constant  system 

x'  =  A  x.  Then  there  exist  M,  a  >  0  such  that  | |X(t)| |  <_  Me"at 

for  t  >  0.  If  x(t;0,x  )  is  a  solution  to  x'  =  B(t)  x  +  b(t), 

then  as  before, 

t 
x(t)  =  X(t)xn  +  /  X(t-s)[(B(s)-A)x(s)  +  b(s)]  ds, 
0    0 

from  which 

l|x(t)||  £  ||x(t)||  ||x  ||  +  r  ||x(t-s)|| 

0     0 

(||B(s)-A||  ||x(s)||  +  ||b(s)||)  ds. 
Then  since  | |b(t) |   :  N 

||x(t)||  eat  <  M  ||x  I  |  +  /  Me"aS||P>Cs)-A||  ||x(s)||d« 

0     0 


t 
0 


+  /  MNe"aS  ds 


<M  ||xl|  +  /  Me-aS||x(s)||  ||B(s)-A||ds 
0     0 

00 

+  /  MNe"aS  ds 


and  thus 

||x(t)||  ea\  L+  /  Me"aS||x(s)||  |  |B(s)-A|  |  ds, 

0 

where  L  is  a  convenient  constant.  Apnlyinq  Gronwall's  Lemma  yields 

t  t 

| | x ( t) | |eaL  <  L  exo  /  M  | |B(s)-A| |  ds, 

0 
so  that 

|x(t)||eau  ■:  L  exo  /  M  |  |B(s)-A|  |  ds  •:  K  <  » 

0 
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where  K  is  another  constant.  Thus  for  t  >_  0,  |.|x(t)||  <_  Ke"a  <_  K, 
which  shows  that  x(t)  and  hence  the  system  x"  =  B(t)  x  +  b(t)  is 
bounded. 

In  addition,  if  y(t;0,y  )  is  any  other  solution  to  this  system, 
||x(t)  -  y(t)||  <  ||x(t)||  +  ||y(t)||  !2Ke-at 
which  shows  that  |x(t)  -  y(t)||-*0  as  t^°°  and  the  system 
x1  =  B(t)  x  +  b(t)  is  also  asymntotically  stable. 

If  the  characteristic  polynomial  of  A  is  not  assumed  stable, 
but  the  system  x'  =  A  x  is  bounded,  then  a  similar  theorem  holds. 

Theorem  4.6.  If  the  system  x(  =  A  x  is  bounded  on  [0,°°),  b(t) 
and  B(t)  are  continuous  on  [0,°°), 


/  |  b(s)  |  ds  <  °°,  and 
0 


f\  |B(s)  -  A| |  ds  <  -  , 
0 

then  the  system  x'  =  B(t)x  +  b(t)  is  bounded. 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  constant 

system  x'  =  A  x.  Then  there  exists  M  >  0  such  that  ||X(t)|   :  M  for 

t  >_  0.  If  x(t;0,x  )  is  any  solution  to  x'  =  B(t)  x  +  b(t),  then  from 

t 
x(t)  =  X(t)xn  +  /  X(t-s)[(B(s)-A)x(s)+b(s)]  ds, 
0    0 


we  have  successively 


t 
x(t)||  ^  M  |  |  x  I  |  +  /  M  ||B(s)-A||  ||x(s)||  ds 
0     0 

t 
+  /  M  |b(s)|  ds 
0 
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lM  ||x  ||  +  /  M  ||B(s)-A||  ||x(s)||  ds 


0 


+  /  M  | |b(s)| |  ds 
0 

t 
:  K  +  /  M| |B(s)-A| |  ||x(s)| |  ds, 
0 


where  K  is  a  suitable  constant.  Applying  Gronwall's  Lemma  yields 

t 

||x(t)||  :  K  exp  /  M  |B(s)-A||  ds 
0 


K  exp  /  M  |B(s)-A|  ds  <  L  <  °° 
0 

where  L  is  another  constant.  Thus  x(t)  and  hence  the  system 

x'  =  B(t)  x  +  b(t)  are  bounded. 

In  general,  if  A(t)  is  not  constant,  Theorems  4.3  through  4.6 

need  no  longer  apply.  Thus  if  a  homogeneous  system  x'  =  A(t)  x 

is  bounded  (stable)  and 

oo 

/  ||b(t)||  dt  <  », 
0 

it  need  not  follow  that  the  non-homogeneous  system  x'  =  A(t)  x  +  b(t) 

is  bounded.  The  same  applies  if  the  homogeneous  system  is  asymptotically 

stable.  A  similar  problem  was  encountered  in  attempting  to  effect  a 

stability  comparison  between  two  non-constant  linear  homogeneous  systems. 

There  it  was  noted  that  additional  hypotheses  were  required  to  insure 

the  stability  of  the  system  x1  =  B(t)  x  when  B(t)  is  asymptotic  to  a 

non-constant  matrix  A(t)  and  the  system  x'  =  A(t)  x  is  stable.  One 

reason  for  this  is  that  stability,  as  defined  in  the  sense  of  Lyaounov, 
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is  quite  delicate  and  often  may  not  be  maintained  even  when  small 
oerturbations  are  introduced.  This  leads  to  the  concepts  of  uniform 
and  restrictive  stability,  which  are  discussed  in  section  V. 
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V.  UNIFORM  AND  RESTRICTIVE  STABILITY 

Since  Lyapunov  stability  does  not  suffice  to  insure  certain 
boundedness  and  comparison  properties,  we  are  led  to  a  consideration 
of  stronger  forms  of  stability.  One  of  these  is  uniform  stability, 
defined  as  follows. 

Definition  5.1 .  Let  x(t;t  ,x  )  be  a  solution  of  the  system 
x'  =  f(t,x),  where  f(t,x)  is  defined  and  continuous  in  the  region 
r=  {(t,x):  t  <  t  <  »,   |x|  <  e,   }.  Then  x(t)  is  said  to  be  uniformly 
stable  if  it  satisfies  (3.2i,ii,  and  iii)  and  in  addition  (5.1) 
given  e  >  0  there  exists  6  =  6(e,f)  >  0  such  that  every  solution 
y(t;tQ,yo)  with  |  |x(*.|J  -  y(t,),J,|  <  $  for  some  t-j  >_  tQ  satisfies 

||x(t)  -  y(t)||  <  e  for  t  >_  tr 

Comnaring  the  definitions,  we  note  that  uniform  stability  is 
similar  to  Lyapunov  stability.  However,  there  is  an  important 
distinction  between  the  two.  A  solution  x(t)  is  stable  in  the  sense 
of  Lyapunov  if  any  other  solution  which  is  initially  close  to  x(t) 
remains  in  an  e-neighborhood  of  x(t)  for  t  >_  t  ;  at  later  times, 
i.e.,  t  >  t  ;  other  solutions  may  come  arbitrarily  close  but  need 
not  remain  there.  If  x(t)  is  uniformly  stable,  once  a  solution  enters 
a  prescribed  small  neighborhood  of  x(t)  for  any  t,  >_  t  ,  it  must  remain 
there  for  all  t  >_  t, ,  Thus,  in  uniform  stability,  6  no  longer 
depends  on  t  . 

The  above  remarks  show  that  a  uniformly  stable  solution  is  stable. 
If  all  solutions  to  a  system  are  uniformly  stable,  we  say  that  the 
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system  is  uniformly  stable.  For  linear  homoqeneous  systems  of  the 

form 

(5.2)  x1  =  A(t)  x, 

if  one  solution  to  the  system  is  uniformly  stable,  then  the  system 

is  uniformly  stable.  Uniform  stability  for  a  system  (5.2)  may  be 

characterized  by  the  behavior  of  the  fundamental  matrix  of  the 

system  as  the  following  theorem  shows. 

Theorem  5.1 .  The  system  x'  =  A(t)  x  is  uniformly  stable  iff 
there  exists  M  >  0  such  that  |X(t)X"  (s)  |  <_  M   t  >_ "s  >_  zn ,  where 
X(t)  is  the  fundamental  matrix  of  the  system. 

Proof:  Assume  that  the  system  x1  =  A(t)  x  is  uniformly  stable, 
so  that  the  trivial  solution  z(t)  =  0  is  uniformly  stable.  Then 
if  e  >  0  there  exists  6  >  0  such  that  if  |x(t,)||  <  6  for  some 
t-,  >_  t  ,  then  |  |x(t)  |  |  <  e  for  t  >_  t, .  Let  y(t;t  ,y  )  be  a  solution 
of  x'  =  A(t)  x.  Then  y(t)  may  be  written 

y(t)  =  X(t)  f](t})   y(t1). 
By  choosing  y(t)  to  be  that  solution  satisfying  y(t,)  =  ■*-  e1 ,  where 

i  A 

e     is  as  in  the  proof  of  Theorem  3.1,  we  have     |y(t, )j|   =  j  <  & 
so  that 

,-1 


|||X(t)X"l(t1)e1||    <  e         tltr 


or  equivalently, 


|  iz"1  (t,t1 )  |  |<  -^1  t  1  t1 

where  Z1  (t,t-. )  is  the  ith  column  of  the  matrix  X(t)X"  (t-,).  Repeating 
the  above  Drocedure  for  i  =  l,...,n  yields 

||x(t)x-1.(t1)||  <  2g£  tLtr 
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Finally,  since  this  procedure  is  valid  for  any  t,  >_  t  and  the 
bound  is  independent  of  t-, ,  we  obtain 

|  |  X ( t )X_1  ( t1 )  |  |  <  ^~   =  M,  t  >_   t1  >_   tQ. 
Conversely,  if  |  |X(t)  X"1  (t-, )  |  |  <  M  for  t  >  tj  >_  t  ,  and 
x(t;t  ,x  )  and  y(t;t  ,y  )  are  any  two  solutions  of  x'  =  A(t)  x, 
then  for  t  >_  t,  >_  t 

||x(t)-y(t)||   =   |  |X(t)X_1  (t^CxCt^-yCt^:)  | 
<_M  \\x{t})   -  yfa)\\ 
so  that  if  e  >  0,  |  |x(t, )  -  y(t, )|  |  <  &  t,  >_t  implies 

|x(t)  -  y(t)|  <  e  t  >  t,  which  shows  that  the  system  is  uniformly 
stable. 

From  Theorem  5.1,  we  obtain  the  following  simple  proposition. 

Proposition  5.1 .  If  the  system  x'  =  A  x,  where  A  is  constant, 
is  stable  on  [0,°°),  then  it  is  uniformly  stable  on  [0,<»). 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  x'  =  A  x.  Then 
there  exists  M  >  0  such  that  ||X(t)||  <_ M  for  t  >  0.  Then  by  (2.13), 

||X(t-s)||  =  |  |X(t)  X_1(s)||  <_M   for  t  >_  s  >_  0 
and  the  system  is  thus  uniformly  stable. 

Recall  that  in  sections  III  and  IV,  several  boundedness  and 
comparison  properties  which  hold  when  A(t)  =  A  is  constant,  need  not 
apply  when  A(t)  is  not  constant.  However,  Proposition  5.1  shows 
that  in  the  case  A(t)  is  constant,  the  systems  considered  are  actually 
uniformly  stable,  rather  than  just  stable.  This  leads  to  the  question 
of  whether  or  not  uniform  stability  of  the  non-constant  system 
x'  =  A(t)  x  is  sufficient  to  insure  these  boundedness  and  comparison 
properties.  The  following  theorems  provide  an  affirmative  answer 
to  this  question. 
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on 


Theorem  5,?.  If  the  system  x1  =  A(t)  x  is  uniformly  stable 

[t  ,°°),  B(t)  is  continuous  on  [t  ,»)  and  ./"I  |  B(t)-A(t)  |  |  dt 

o 

then  the  system  x'  =  B(t)  x  is  uniformly  stable. 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  system 

x'  =  A(t)  x.  Then  |  |X(t)  X"1  (s )  |  |  <  M  for  t  >_  s  >_  t  .  Let  t   >  0 

and  fix  s  >  t  .  If  x(t;t  ,x  )  and  v(t;t  ,y  )  are  two  solutions  of 
—  o         oo     -oo 

the  system  x1  =  B(t)x,  then  since  B(t)  =  A(t)  +  [B(t)  -  A(t)], 
we  have  for  t  >_  s, 

x(t)-y(t)  =  X(t)X-1(s)[x(s)-y(s)] 

+  /  X(t)X-1(T)[B(x)-A(x)][x(T)-y(T)]  dx 


and  thus 


x(t)-y(t)||  <_  ||X(t)X"1(s)||  ||x(s)-y(s)|| 

t  1 

+  /  II  X(t)X"'(T)||  | |B(t)-A(x) I I I |x(x)- 


y(x)| |dx 


M  ||x(s)-y(s) 


o 
Applying  Gronwall's  Lemma  yields 


t 
+   /  M  ||B(t)-A(t)||  I |xCx)-yCx) | |  dx. 
t 


from  which 


t 
x(t)-y(t)||  :  M  ||x(s)-y(s)||  exn  /  M||B(x)-A(x) 

'°        d 
dx. 


x(t)-y(t)||  :  M  ||x(s)-y(s)||  exD  /  M  ||B(x)- 

A(x)||  dx. 
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Therefore, 

||x(t)-y(t)||  <_L  ||x(s)-y(s)|| 
where  L  is  a  convenient  constant.  The  above  relation  is  valid 
for  t  ^  s  >_  t  ,  so  that  if  |x(s)-y(s)|  <  r-    s  >_  t  ,  then 

| x ( t )-y ( t ) |  <  e  for  t  ^s  which  shows  that  the  system  x'  =  B(t)  x 
is  uniformly  stable. 

Corollary  1 .  Under  the  same  hypotheses  as  Theorem  5.2,  the 
system  x'  =  B(t)  x  +  b(t),  where  b(t)  is  continuous  on  [t  ,»)  is 
uniformly  stable  on  [t  ,°°). 

Proof:  The  proof  follows  from  the  fact  that  under  the  hyootheses, 
the  system  x'  =  B(t)  x  is  uniformly  stable  and  that  if  x(t)  and 
y(t)  are  two  solutions  to  the  non-homogeneous  system  x*  =  B(t)  x  +  b(t), 
then  the  difference  x(t)  -  y(t)  is  a  solution  to  the  non-homogeneous 
system  x'  =  B(t)  x. 

When  the  homogeneous  system  is  uniformly  stable,  then 
boundedness  of  the  non-homogeneous  system  may  be  determined  as 
for  the  constant  case  considered  earlier  in  Section  IV. 

Theorem  5.3.  If  the  system  x1  =  A(t)  x  is  uniformly  stable 

on  [t  ,°°),  b(t)  is  continuous  on  [t  ,°°)  and   /1°|b(s)  |  |  ds  <  °°, 

o  °       0 

then  the  system  x'  =  A(t)  x  +  b(t)  is  bounded  on  [t  ,°°). 

Proof:  Let  x(t;t  ,x  )  be  a  solution  to  x'  =  A(t)  x  +  b(t). 

Then 

x(t)  =  X(t)x  +  /  X(t)X"' (s)b(s)  dss 

Jo 

where  X(t)  is  the  fundamental  matrix  of  the  homogeneous  svstem 

x'  =  A(t)  x.  Since  |  |X(t)  X"1  (s)  |  |  <_  M  for  t  >_  s  >  t  , 
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t 

x(t)  I  I  <_  M  I  |x  |  |  +  /  M  I  |b(s)|  |  ds 
0     t 


and  thus 


|x(t)  |  <_  M  |  | x  |  +   /  M  |b(s)  |  |  ds  <_  L<  oo 

to 
where  L  is  a  convenient  constant.  Thus  the  system  x'  =  A(t)  x  +  b(t) 

is  bounded. 

It  is  also  possible  to  determine  boundedness  of  a  non-homoqeneous 
system  through  comparison  with  a  uniformly  stable  system. 

Proposition  5.2.  If  the  svstem  x'  =  A(t)  x  is  uniformly  stable 

oo 

on  [t  ,~),  B(t)  and  b(t)  are  continuous  on  [t  ,°°),   /||b(t)|  dt  <  °°, 
o  o     t 

o 

oo 

and   /||B(t)  -  A(t) | |  dt  <  »,  then  the  system  x*  =  B(t)  x  +  b(t) 

*o 

is  bounded  on  [t  ,°°) . 
o 

Proof:  By  Theorem  5.2,  the  hypotheses  imply  that  the  homo- 
geneous system  x'  =  B(t)  x  is  uniformly  stable  on  [t  ,°°) .  The 
result  then  follows  from  Theorem  5.3. 

Thus  certain  boundedness  and  comparison  properties  established 
for  constant  systems  may  be  generalized  to  non-constant  systems 
which  are  uniformly  stable. 

We  now  introduce  the  concept  of  restrictive  stability. 

Definition  5.2.  The  linear  homogeneous  system  x'  =  A(t)  x 
is  said  to  be  restrictively  stable  if  it  is  stable  and  its  adjoint 
system  x'  =  -A(t)  x  is  stable. 

We  note  that  restrictive  stability  applies  to  a  system,  rather 
than  a  particular  solution  and  is  defined  for  linear  homogeneous 
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systems  only.  The  following  theorem  characterizes  restrictive 
stability  in  terms  of  the  behavior  of  the  fundamental  matrix  of  the 
system. 

Theorem  5.4.  The  system  x1  =  A(t)  x  is  restrictively  stable 
on  [t  ,-)  iff  there  exists  M  >  0  such  that  ||X(t)||  ^M  and 
I  I  X"  (t)||<_M  for  t  >_  t  ,  where  X(t)  is  the  fundamental  matrix 
of  the  system  x1  =  A(t)x. 

Proof:  The  result  follows  by  Theorem  3.1,  noting  that  [X"  (t)] 
is  the  fundamental  matrix  of  the  adjoint  system  x'  =  -A(t)  x. 

From  Theorem  5.4,  we  see  that  a  restrictively  stable  system 
x'  =  A(t)  x  is  uniformly  stable  since 

I  [  X(t)X"1  (s)  t  I  1  ||X(t)||  I  I  X"1  Cs)  1 1  <_M2     ».t>'t0 
but  not  necessarily  conversely.  Also,  since  -[-A(t)  ]  =  A(t), 
it  follows  that  if  x'  =  A(t)  x  is  restrictively  stable,  then  its 
adjoint  svstem  is  restrictively  stable.  The  following  propositions 
give  conditions  which  suffice  for  a  system  to  be  restrictively 

stable. 

00 

Proposition  5.3.  If   /  ||A(t)|  dt  <  <*>,  and  A(t)  is  continuous 
on  [t  ,«=),  then  the  system  x1  =  A(t)  x  is  restrictively  stable  on 

[t0.-). 

Proof:  Proposition  3.2  shows  that  under  the  hypotheses,  the 
system  x'  =  A(t)  x  is  stable.  Since  ||A(t)||  =  ||  -A(t)T||,  it 
follows  that  the  adjoint  system  is  stable. 

Proposition  5.4.  If  the  system  x1  =  A(t)  x  is  stable  on  [t  ,°°) 

and  A(t)  is  continuous  on  [t  ,«>),  and  lim  inf  /  tr  A(s)  ds  >  -°°, 

t+°°      t 
0 

then  it  is  restrictively  stable. 
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Proof:  Let  X(t)  be  the  fundamental  matrix  of  x'  =  A(t)  x.  Then 
X(t)  is  bounded  for  t  >_  t  since  the  system  is  stable.  By  Lemma  3.3, 

lim  inf   /  tr  A(s)  ds  >  -°°  implies  that  X"  (t)  is  bounded  for 

t+°°      t 
o 

t  >_  t  .  Hence  the  system  x1  =  A(t)  x  is  restrictively  stable  by 
Theorem  5.4. 

Proposition  5.5.  If  the  adjoint  system  x'  =  -A(t)  x  is  stable 

on  [t  ,°°),  A(t)  is  continuous  on  [t  ,°°),  and  lim  sud  /  tr  A(s)  ds  <  °°, 

0  °        t-*»      t 

o 

then  the  system  x'  =  A(t)  x  is  restrictively  stable  on  [t  ,°°). 

o 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  x '  =  A(t)  x.  Then 
since  the  adjoint  system  is  stable,  X"  (t)  is  bounded  for  t. >,  t  . 
From  the  relation  X(t)  =  det  X(t)  adj  X-1(t),  we  obtain  by  (2.8) 

1  1 

X(t)  =  exp  /  tr  A(s)  ds  adj  X~'(t) 

t 
o 

and  since  X"  (t)  is  bounded,  adj  X"  (t)  is  bounded  and  thus 

t 

||X(t)||  <_M  exn  r   tr  A[s)  ds 

t 
:  M  lim  sup  /  tr  A(s)  ds  <   L  <  °° 

t+oo  t 

0 

where  M  and  L  are  suitable  constants.  Hence  X(t)  is  bounded  and 
the  system  x1  =  A(t)  x  is  restrictively  stable  by  Theorem  5.4. 

An  interesting  result  concerning  the  well  known  second  order 
Hill ' s  Equation 

(5.3)  u"  +  a(t)  u  =  0   on  [0,»), 

where  a(t)  is  real  valued  and  continuous  on  [0,°°)  and  u(t)  is 
a  scalar  valued  function,  may  be  obtained  through  the  use  of 
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stability  theory.  The  vector  analogue  of  (5.3)  is  the  two-dimensional 

system 

(5.4)  y'  =  z 

z'  =  -a(t)  y  J  . 

which  is  in  the  form  (5.2)  with  x  =  (y)     and  A(t)  =       |. 

VZ/  \-a(t)  0/ 

We  note  that,  since  tr  A(t)  e  0,  if  the  system  (5.4)  is  stable, 

then  it  will  be  restrictively  stable,  and  thus  uniformly  stable. 

00 

It  is  now  possible  to  qive  the  followinq  result:  If   /|a(t)|  dt  <  », 

0 

then  there  exist  unbounded  solutions  of  (5.3).  This  may  be 

proved  using  stability  theory  as  follows.  Assume  that  all  solutions 

to  (5.3)  are  bounded.  Then  (5.4)  is  a  stable  system  and  hence 

uniformly  stable  by  the  above  remarks.  Theorem  5.2  then  asserts 

that  if 

/O    1 

B(t)  = 

\0    OJ 

then  the  system  x'  =  B(t)  x  is  uniformly  stable  and  hence  bounded 

oo  oo 

since   /  ||B(t)  -  A(t)||  dt  =   /  |a(t)|  dt  <  oo.  But  this  is  a 
0  0 

contradiction  since 

*'  ;o 

is  a  solution  to  x'  =  B(t)x.  Thus  all  solutions  to  (5.3)  are  bounded. 

Restrictive  stability  may  also  be  characterized  in  terms  of 
the  reducibility  of  the  system.  Toward  this  end,  we  have  the 
following  definition. 
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Def in : Lion  5.3.  The  system  x*  =  A(t)  x  is  said  to  be  reducible 
if  there  exists  a  matrix  L(t)  whose  elements  are  differentiate 
functions  which  is  bounded  together  with  L"  (t)  on  [t  ,»)  such 
that 

(5.5)  L'^t)  A(t)  L(t)  -  L_1(t)  L'(t)  =  B 
where  B  is  a  constant  matrix.  If  B  =  0,  the  zero  matrix,  then 
we  say  that  the  system  is  reducible  to  zero. 

The  relation  (5.5)  is  often  expressed  by  sayinq  that  A(t)  is 
kinematically  similar  to  the  constant  matrix  B.  If  both  L(t)  and 
A(t)  are  constant,  then  (5.5)  becomes 

(5.6)  L"1  A  L  =  B 

which  is  the  ordinary  similarity  relation. 

If  a  linear  homogeneous  system  is  reducible,  then  the  trans- 
formation x  =  L(t)  y  reduces  the  system  xf  =  A(t)  x  to  the  constant 
system  y'  =  By,  while  if  x1  =  A(t)  x  is  reducible  to  zero,  then 
this  transformation  yields  the  system  y'  =0. 

The  following  theorem  characterizes  restrictive  stability  in 
terms  of  the  reducibility  of  the  system. 

Theorem  5.5.  The  system  x'  =  A(t)  x  is  restrictively  stable  on 

[t  ,°°)  iff  it  is  reducible  to  zero  on  [t  ,°°). 
L  o  o 

Proof:  If  x1  =  A(t)  x  is  restrictively  stable,  then  there 
exists  M  >  0  such  that  |  |X(t)  1 1  <_  M  and  |  |X_1  (t)  |  |  <  M  f or  t  >_  t  , 
where  X(t)  is  the  fundamental  matrix  of  the  system.  Then  since 
X(t)  satisfies  X(t)  =  A(t)  X,  it  follows  that 

X~](t)  A(t)  X(t)  -  X_1(t)  X'(t)  e  0 
and  that  the  system  x'  =  A(t)  x  is  reducible  to  zero  by  taking 
L(t)  =  X(t). 
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Conversely,  assume  that  (5.5)  is  satisfied  for  B  =  0.  Then 
from  (5.5),  L(t)  satisfies 

L'  =  A(t)  L 
and  since  L(t)  is  non-singular,  it  follows  that  L(t)  is  a  basis 
matrix  for  the  system  x'  =  A(t)  x.  Hence  X(t)  =  L(t)  W,  where  X(t) 
is  the  fundamental  matrix  and  W  is  a  constant  non-singular  matrix. 
Thus  ||X(t)||  <_  ||L(t)||  !|W||  and  |  | X"1  (t)  |  |  <  ||Lfc,(t)||  1 ( W"1 1 1  which 
shows  that  both  X(t)  and  X"  (t)  are  bounded.  Thus  the  system 
x1  =  A(t)  x  is  restrictively  stable  by  Theorem  5.4. 

Uniform  stability  is  also  related  to  reducibility  as  the 
following  theorem  illustrates. 

Theorem  5.6.  If  the  system  x'  =  A(t)  x  is  stable  and  reducible 
on  [0,°°),  then  it  is  uniformly  stable  on  [0,°°). 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  x'  =  A(t)  x  and 
let  L(t)  be  as  in  Definition  5.3.  Then  the  transformation 
y(t)  =  L"  (t)  x(t)  reduces  x1  =  A(t)  x  to  the  system  y*  =  B  y, 
where  B  is  constant.  Also,  since  L~  (t)  X(t)  is  non-singular  and 


satisfies 


U"1  X)'  =  B  (L_1  X), 


it  follows  that  L~  (t)  X(t)  is  a  basis  matrix  for  y'  =  B  y.  Then 
Y(t)  =  L_1(t)  X(t)  W,  where  Y(t)  is  the  fundamental  matrix  of 
y1  =  B  y  and  W  is  a  constant  non-singular  matrix.  Since  x'  =  A(t)  x 
is  stable,  X(t)  is  bounded  and  by  definition,  both  L(t)  and  L"  (t) 


are  bounded.  Thus  there  exists  M  >  0  such  that  |X(t)|   :  M, 
||L(t)||  <_M,  and  ||L_1(t)||  <_M,  so  that||Y(t)||  <_  M2  ||W||,  which 
shows  that  Y(t)  is  bounded  and  thus  that  the  system  y'  =  B  y  is  stable. 
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Then  by  Proposition  5.1,  y'  =  B  y  is  uniformly  stable  and  thus  there 
exists  N  >  0  such  that  |  |Y(t)  Y"1  (s)  |  |  <_  N  for  t  >_  s  ^  0.  Then 
since 

X(t)  =  L(t)  Y(t)  W"1 
we  have  for  t  >_  s  >_  0 

||X(t)X-1(s)||  <  ||L(t)||  ||Y(t)Y-1(s)||  ||L"J(t)|| 
<_  M2  N 
and  thus  x'  =  A(t)  x  is  uniformly  stable  by  Theorem  5.1, 

The  above  remarks  relate  the  connection  between  uniform  and 
restrictive  stability  in  a  system  and  the  reducibility  of  that 
system.  Also,  we  have  seen  that  uniformly  stable  systems  and 
restrictively  stable  systems  satisfy  certain  boundedness  and 
comDarison  properties  that  need  not  necessarily  hold  for  stable 
systems.  In  Section  VI,  we  shall  show  that  a  similar  behavior 
is  exhibited  when  the  systems  considered  are  non-1 i near. 
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VI.  NONLINEAR  SYSTEMS 

The  systems  considered  here  are  non-linear  mixed  systems  of 
the  form 

(6.1)  x'  =  A(t)  x  +  f(t,x) 

where  A(t)  is  continuous  on  [0,°°)  and  f(t,x)  satisfies  certain 
aonroDriate  conditions.  Such  systems  are  extremely  important 
since  the  general  first  order  normal  system 

(6.2)  x'  =  g(t,x) 

can  often  be  put  in  the  form  (6.1)  by  taking  f(t,x)  =  q(t,x)  -  A(t) 
where  A(t)  depends  on  the  form  of  g(t,x).  Conditions  on  f(t,x) 
are  thus  sought  which  insure  stability  in  the  non-linear  system 
(6.1)  when  it  is  known  that  the  linear  portion 

(6.3)  x'  =  A(t)  x 

is  stable.  When  A(t)  =  A  is  constant,  then  (6.1)  becomes 

(6.4)  x'  =  A  x  +  f(t,x) 
and  stability  is  easier  to  determine. 

Consider  the  following  conditions  on  f(t,x). 
(6.5i)  f(t,x)  is  continuous  in  the  region 

r  =  {(t,x):  0  <  t  <  co,  ||x|  |  <  U 
(6.5ii)  |f(t,x) | |  =  o  ( |  |x|  | )  uniformly  in  t  as 

|x||+0,  i.e.,  there  exists  a  function  u(r)  >  0 
such  that  u(r)+0  as  r+0  and 
l|f(t,x)||  <  ||x||  u(||x||). 
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Condition  (6.5i)  assures  existence  of  solutions  for  small 
|x|  but  not  uniqueness.  Note  also  that  (6.5ii)  imDlies  that 
z(t)  =  0  is  a  solution  to  (6.4). 

When  f(t,x)  satisfies  (6.5i,ii),  we  might  heuristically 
reason  that  since  f(t,x)  is  small  when  |x|  is  near  zero,  the 
system  (6.4)  miqht  be  expected  to  behave  in  a  manner  similar 
to  x1  =  A  x.  Such  reasoning  is  correct  as  the  following  theorem 
illustrates. 

Theorem  6.1 .  If  the  characteristic  Dolynomial  of  A  is  stable 

and  f(t,x)  satisfies  (6.5i,  ii)  then  the  solution  z(t)  =  0  of 

the  system  x'  =  A  x  +  f(t,x)  is  asymptotically  stable. 

Proof:  We  first  show  that  solutions  x(t;0,x  )  exist  and  are 

o 

defined  on  [0,°°)  when  |x  |  is  small.  Let  X(t)  be  the  fundamental 

matrix  of  the  constant  system  x1  =  A  x.  Then  X(t)  satisfies  for  t  >_  0 

|X(t)|  <  M  e"a  for  suitable  M  and  a.  The  solution  x(t;0,x  )  of 
—  o 

x1  =  A  x  +  f(t,x)  satisfies  (where  it  exists) 

t 
x(t)  =  X(t)  x  +  /  X(t-s)  f(s„x(s))  ds, 
0    0 

so  that 

||x(t)||eat  <..M  ||xl|  +  /  MeaS||f(s,x(s))||  ds. 

0 

The  above  inequality  is  valid  for  t  [0,t-,)  where  |x(t)||  <  £for 

t  [0,^),  assuming  ||xo||  <  &     Then  from  (6.5ii),  given  E  >  0,  there 

exists  6  >  0  such  that  for  t  >  0, 

l|f(t,x)||  <  e||x|| 
whenever  |x|  <  6. 


52 


If  we  take  | |x  |  <  6,  then  by  continuity  of  x(t),  there  exists 
t2  >  0  such  that  |  |x(t)  |  |  <  &   for  0  <_  t  <  t„.     Therefore  for 
0  <_  t  <  T,  where  T  =  min  (t-,  ,t«)  •  we  have 

||x(t)||  eat  <  M  ||x  ||  +  Me  /  eaS||x(s)M  ds. 

0        0 

Applyinq  Gronwall's  Lemma  to  the  above  yields 

'  t  t 

||x(t)| |eaL  <  M| |x J  exp  /  Me  ds 

0       0 
and  thus  for  0  <_  t  <  T 

I |x(t)| |  1  M  | |x  | |  exp(Me-a)t. 
Choosing  e  <  £■  yields 

||x(t)||  V'M  ||x0||   0  <_t  <.T 
so  that  if  ||xo|!  <  -^  , 

|  |x(t)|  |  <  |  <  6   0  <_  t  <_  T. 

Since  f(t,x)  is  defined  and  continuous  for  t  >_  0  and  |x|  <  5, 

the  above  relation  imDlies  that  the  solution  x(t;0,x  )  can  be 

o 


extended  interval  by  interval  Dreservinq  the  above  bound.  Thus, 

<0)  with  MX0II  <  a| 


any  solution  x(t;0,x  )  with  |x  |  <  m  exists  on  [0,°°)  and  satisfies 


|x(t)  I  <  j    for  t  >_  0.  Since  6  can  be  taken  arbitrarily  small , 
this  implies  that  the  solution  z(t)  =   0  of  x1  =  A  x  +  f(t,x) 
is  stable.  Furthermore, 

I  1  x(t)  |  |  <_  M  |  |x  |  |  exo  (Me-a)t 
where  e  <  S-  implies  that  |x(t)||->0  as  t+°°  and  thus  the  trivial 
solution  is  asymptotically  stable. 
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Coddington  and  Levinson  [3]  show  that  the  hypotheses  of  Theorem  6.1 
may  be  weakened  somewhat.  Instead  of  requiring  that  f(t,x)  satisfy 
(6.5ii)  it  is  sufficient  to  assume  instead  that  for  some  K  >  0, 

|  |f(t,x)|  |  <_  K  |  |x|  |  t  >_  0 
for  all  small   |x||,  and  that  qiven  e  >  0  there  exist  6,T  such  that 

IK(t,x)||  ie  ||x|| 
for  |x|  <  6  and  t  >_  T. 

Theorem  6.1  may  then  be  applied  to  a  system  of  the  form 
(6.6)  x'  =  A  x  +  B(t)  x  +  g(t,x) 

where  B(t)  -»-0  as  t-*»  and   ||g(t,x)||  =  o(||x||). 

Other  conditions  on  f(t,x)  which  suffice  for  theorems  similar 
to  Theorem  6.1  may  be  found  in  [1]  and  [6]. 

When  A(t)  is  not  constant  in  (6.1),  then  Theorem  6.1  need  no 
longer  necessarily  hold.  However,  if  the  linear  portion  xf  =  A(t)  x 
of  (6.1)  is  uniformly  stable,  then  we  have  the  following  theorem. 

Theorem  6.2.  If  the  linear  system  x1  =  A(t)  x  is  uniformly  stable 
and  f(t,x)  satisfies 
(6.7i)  f(t,x)  defined  and  continuous  for  t  >_  0, 

|xj  |  <£  and 
(6.7ii)  ||f(t,x)||  <_a(t)  ||x||,  where  a(t)  >  0  and 

oo 

/  a(t)  dt  <  oof 
0 

then  the  solution  z(t)  e  0  of  x'  =  A(t)  x  +  f(t,x)  is  uniformly 
stable  on  [0,°°). 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  linear  system 
x'  =  A(t)  x.  Then  there  exists  M  >  0  such  that  |  |X(t)  X_1(s)||  ^M 
for  t  >_  s  >_  0.  As  in  the  proof  of  Theorem  6.1,  it  can  be  shown  that 
any  solution  x(t;0,xo)  of  x'  =  A(t)  x  +  f(t,x)  exists  and  is  defined 
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on  [0,«>).  Then  assuming  this,  the  solution  x(t;0,x  )  can  be  written 

1         *      1 

X(t)  =  X(t)X"'(s)x(s)  +  /   X(t)X~'(x)f(x,x(x))  dx, 

0 

where  s  >_  0  is  fixed,  so  that 

t 
||x(t)||  <  M||x(s)||  +  /  M  a(T)| |x(t)| |dx. 

0 

Applvinq  Gronwall's  Lemma  yields 

t 

| |x(t)| [<  M| |x(s)| |  exn  /  M  a(x)  dx 

0 

and  thus 

oo 

||x(t)||  :  M||x(s)||  exp  /  M  a(x)  dx  -:  K||x(s)|| 

0 

where  K  is  a  convenient  constant.  The  above  relation  is  valid  for  any 

s  >_  0  so  that  given  e  >  0,  ||x(s)||  <  i-  for  some  s  >_  0  implies 

|x(t)||  <  e  for  all  t  >_  s  so  that  the  trivial  solution  z(t)  =  0  of 

x*  =  A(t)  x  +  f(t,x)  is  uniformly  stable. 

Corollary  1 .  Under  the  same  hypothesis  as  Theorem  6.2,  if  the 
linear  system  x'  =  A(t)  x  is  in  addition  asymptotically  stable,  then 
so  is  the  trivial  solution  z(t)  =  0  of  x'  =  A(t)  x  +  f(t,x). 

Proof:  The  hypotheses  imply  that  z(t)  =  0  is  a  uniformly  stable 

solution  of  x'  =  A(t)  x  +  f(t,x).  If  the  linear  system  is  in  addition 

asymptotically  stable,  then  ||X(t)||  ^0  as  t->~  where  X(t)  is  the 

fundamental  matrix  of  the  linear  svstem.  Therefore,  if  x(t;0,x  )  is  anv 

o 

solution  to  x'  =  A(t)  x  +  f(t,x),  then  aiven  e>0  there  exists  T  such 

that  ||X(t)||  <  |-  for  t  >  T.  Then  from 
xo 

1      1 
x(t)  =  X(t)xA  +  /  X(t)X  '(s)f(s,x(s))ds, 

0   0 


we  have 


x(t)||  i  ||X(t)||  ||x  l|  +  /  M  a(s)  ||x(s)||  ds, 

0    0 
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where  |  |X(t)  X"1  (s)  |  |  <_  M  for  t  >_  s  >_  0.  Then  for  t  >_  T  we  obtain 

t 
||x(t)||  <  e  +  M   /  ats)  ||x(s)||  ds 
0 
and  applying  Gronwall's  Lemma, 

t 
| |x(t)| |  :  e  exp   ;  M  a(s)  ds, 
0 

oo 

so  that  ||x(t)||  :  £  exp  /  M  a(s)  ds  <  eK 

0 

where  K  is  a  suitable  constant.  Since  e  >0  is  arbitrary  and  K  is 

independent  of  e  and  T,  we  conclude  that  |  |x(t)||->0  as  t-*»  and 

that  the  trivial  solution  z(t)  =  0  of  x'  =  A(t)  x  +  f(t,x)  is 

uniformly  and  asymptotically  stable. 

From  the  proof  of  Theorem  6.1,  we  note  that  when  the  linear 
system  x'  =  A(t)  x  is  uniformly  stable  and  f(t,x)  satisfies 
(6.7i,ii),  then  all  solutions  of  (6.1)  remain  bounded.  However, 
in  the  event  the  linear  system  is  restrictively  stable,  then  the 
solutions  of  (6.1)  are  actually  bounded  by  those  of  the  linear 
system  as  the  following  theorem  illustrates. 

Theorem  6.3.  If  the  system  x'  =  A(t)  x  is  restrictively  stable 
and  f(t,x)  satisfies  (6.7i,ii),  then  ewery   solution  x(t;0,x  ) 
of  x1  =  A(t)  x  +  f(t,x)  satisfies  x(t)  =  0(||X(t)||),  i.e., 
| |x(t)| |  <_  K  | |X(t)| |  | |x  |  |  where  K  >  0  is  constant. 

Proof:  Let  X(t)  be  the  fundamental  matrix  of  the  linear  system 
x'  =  A(t)  x.  Then  there  exists  M  >  0  such  that  ||X(t)||  _^M  and 
|  |X-1(t)|  |  <_  M  for  t  >  0.  Let  x(t;0,x  )  be  a  solution  of 
x'  =  A(t)  x  +  f(t,x).  Then  from 

x(t)  =  X(t)x  +   /  X(t)  X^'(s)  f(s,x(s))  ds, 

0    0 
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we  obtain 


x(t) 

x(t) 

1  I |x  I  I  +  M  /  a(s)  |  |x(s)||  ds 
0      0 


||x  ||  +  r  /  a(s)  -r££U   ds 
0      0 


Applying  Gronwall's  Lemma  yields 


x(s) 

X(s) 

x(t) 

x(t) 

| |x  I  I  exn  /   M  a(s)  ds 
0      0 


I Ix  II  exn  /  M  a(s)  ds  =  K|  |x  |  | 
0      0         -    o  » 

where  K  is  a  convenient  constant.  Thus  1 1  x(t))ll  IK  Hxll  I  lx(t)  1 1 
which  completes  the  proof. 

Theorem  6.3  shows  that  when  the  linear  system  is  restrictively 
stable,  the  contribution  of  the  non-linear  term  f(t,x)  to  the  solutions 
of  (6.1)  is  small.  Golomb  [4]  has  qiven  conditions  under  which  bounded- 
ness  and  stability  may  be  determined  in  systems  of  the  form  in  which 
f(t,x)  satisfies 


(6.81) 


(6.8ii) 


f(t,x)  is  continuous  in  x  a.e.  for  t  >_  0 
and  is  inteprable  in  t  for  fixed  x  and 

|  |f(t,x)|  |  <_a(t)  u(|  |x|  |)  for  t  >_0,  where 
a(t)  is  positive  and  Intearable  and  u(r) 
defined  for  r>  0  is  a  positive  non-decreasina 
continuous  function  which  is  submultiplica- 
tive,  i.e. , 
u(r  s)  :  u(r)  u(s) , 
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without  requiring  uniform  stability  in  the  linear  system,  essentially 
qeneralizinq  Theorems  6.2  and  6.3.  Other  evaluations  of  solutions  to 
(6.1)  may  be  found  in  the  same  paper. 
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VII.  CONCLUSIONS 

We  have  considered  stability  in  the  sense  of  Lyapunov  for 
systems  of  ordinary  differential  equations  with  particular  emohasis 
on  linear  systems.  The  theory  for  these  systems  is  quite  advanced 
today  and  in  Sections  III,  IV,  and  V  we  have  considered  several  of 
the  more  important  results.  Though  by  no  means  exhaustive,  much  of 
the  material  in  these  sections  comes  from  varied  sources  and  it  is 
felt  that  the  presentation  there  offers  a  more  complete  treatment 
than  is  usually  found  elsewhere.  In  addition,  a  new  set  of  sufficient 
conditions  for  a  stability  comparison  theorem  have  been  formulated 
which  to  the  best  of  the  author's  knowledge  have  not  appeared  else- 
where. Though  the  treatment  in  Section  VI  is  somewhat  brief,  it  is 
intended  for  the  most  part  to  complement  the  rest  of  the  paper  by 
illustrating  some  of  the  more  useful  properties  exhibited  by  uniformly 
and  restrictively  stable  systems.  The  theory  of  stability  for  non- 
linear non-autonomous  systems  is  still  quite  incomplete  and  is  today 
an  area  of  extensive  active  research. 

Some  other  aspects  of  stability  theory  not  considered  here  are 
the  Poincare-Bendixson  theory  for  autonomous  systems,  the  second 
method  of  Lyaounov,  and  various  analytical  approaches  such  as  the 
method  of  Lindstedt,  Krylov  and  Bogolyubov,  and  of  Van  Der  Pol. 
As  mentioned  earlier,  the  approach  to  the  problem  of  stability  of 
a  system  of  differential  equations  varies  with  both  the  form  and 
order  of  the  system.  This  paper  has  considered  but  one  of  these 
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approaches,  one  which  is  particularly  useful  in  the  treatment  of 
linear  systems.  For  a  treatment  of  some  of  the  other  areas  in  the 
theory  of  stability  of  ordinary  differential  equations,  we  refer  the 
reader  to  several  of  the  excellent  texts  on  stability,  including 
[2],  [5]  and  [7]. 
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